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On the basis of the phenomenological relativistic harmonic models for quarks we have obtained the 
masses of P wave mesons. The full Hamiltonian used in the investigation has Lorentz scalar + vector 
confinement potential, along with one gluon exchange potential (OGEP) and the instanton-induced 
quark-antiquark interaction (III). A good agreement is obtained with the experimental masses. The 
respective role of III and OGEP for the determination of the meson masses is discussed. 



PACS numbers: 12.39.Ki, 12.39.Pn, 12.39.-x 



I. INTRODUCTION 



There is a wealth of experimental data in hadron spec- 
troscopy that would constitute a good testing ground 
for nonpcrturbative quantum chromodynamics (QCD). 
Since the exact form of confinement from QCD is not 
known, one has to go for phenomenological models. 
The phenomenological models are either non-relativistic 
quark models (NRQM) 0, S, 0, i, [1 [1, Q with a suitably 
chosen potential, or relativistic models where the inter- 
action is treated perturbatively. NRQM have proved to 
be quite successful in describing the hadronic properties. 
The Hamiltonian of these quark models usually contains 
three main ingredients: the kinetic energy, the confine- 
ment potential and an hyperfine interaction term, which 
has often been taken as an effective one-gluon-exchange 
potential (OGEP) Q. Other types of hyperfine in- 
teraction have been introduced in the literature; from 
the non-relativistic reduction of the t'Hooft interaction 
S EH' termed as Instanton-induced Interaction (III), 
which has already been successf ully applied in several 
studies of the hadron spectra @, El Gj, El El EE EH . 
The main achievement of the III in hadron spectroscopy 
is the resolution of the J7a(1) problem, which leads to 
a good prediction of the masses of rj and rj mesons. 
The Goldstone-Boson-Exchange interaction introduced 
by Glozman and Riska [17| furnishes another example 
of hyperfine interaction; it allows a good description 
of the baryon spectrum and gives a correct order- 
ing for the positive and negative parity states. The 
model of Glozman and Riska has however the major 
caveat to apply only to baryons and is thus not able 
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to give an unified description of the spectrum of hadrons. 

The successes of the NRQM in describing the spec- 
trum of light hadrons are somehow paradoxical, as light 
quarks should in principle not obey a non-relativistic dy- 
namics. This paradox has been avoided in many works 
based on the constituent quark model by using for the 
kinetic energy term of the hamiltonian a semi-relativistic 
or relativistic expression (see, for example, [HI, EH EH) • 

In literature there are models which have tried to ex- 
plain hadron spectroscopy only with OGEP [l], 0] and 
some models with only III 0, EH> ignoring completely 
the OGEP. We feel that it may be an exaggeration to 
eliminate OGEP completely for light quarks. The OGEP 
still has to be present but with a smaller strength con- 
sistent with the asymptotic freedom, since the III has to 
vanish for heavy quarks. 

In our present work, we have made use of the rela- 
tivistic harmonic model (RHM) [20]. The RHM com- 
bined with OGEP has already been used to calculate light 
hadrons masses, baryons magnetic moments, leptonic de- 
cay widths and N-N scattering phase shifts f2lL I22I [2^]. 
Previously, by employing the RHM [H, |24{ and the 
NRQM 25[ along with III, the ground state masses of 
pseudoscalar and vector mesons were investigated. In 
both the cases the results showed that the inclusion of III 
diminished the relative importance of OGEP for the hy- 
perfine splitting. One of the aims of the investigation was 
also to test whether quark gluon coupling constant(a s ) 
can be treated as a perturbative effect and to understand 
the role played by the III in meson spectra [24|. 

In view of the apparent success of RHM [24[ in the 
description of S"-wave spectra of mesons, we feel it is 
worthwhile to apply it to the case of orbitally excited 
states. In our present work the full hamiltonian used 
has a Lorentz scalar plus vector confinement potential, 
along with central and non-central (spin-orbit and ten- 
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sor) terms of OGEP and III. In addition, III also has an- 
tisymmetric spin-orbit term proportional to L ■ A where 
A is defined in terms of the Pauli matrices as ^(c?i — o^)- 
The full discussion of the hamiltonian is given in section 
II. We also discuss the parameters involved in our model 
in section IIIII The results of the calculation are pre- 
sented in section IIVI and the conclusions are given in 
section [V] 

In our present work, we have computed the masses of 
the singlet and triplet light P wave mesons by including 
the III as a short-range nonperturbative gluon effect in 
addition to the perturbative conventional OGEP derived 
from QCD. This will allow much better understanding 
of the P-wave meson spectroscopy, where some of the qq 
quark model assignments of the known mesons are still 
controversial. We hope it will also allow us a better un- 
derstanding of the production properties of the P-wave 
mesons. In literature there are numerous attempts to 
understand the P-wave meson spectroscopy. The refer- 
ence can be found in the review (26|. One of the aims 
of this study is to determine explicitly the role played by 
instantons in meson spectra, when used in the framework 
of the RHM and to compare the effects of III to that of 
OGEP. 
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Here U is a momentum and eigenvalue E dependent 
transformation operator. With this transformation, the 
upper component $ satisfies the equation 
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<f> = (E-M)<Z> 



(7) 



which is like the three dimensional harmonic oscillator 
equation with an energy dependent parameter fl^ 



n 2 n = A{E n + M) 2 
The eigenvalue of ([7]) is thus given by 
E 2 n = M 2 + (2n + l) n 2 n 



(8) 



(9) 



Note that Eq.® can also be derived by eliminating 
the lower component of the wave function, using a Foldy- 
Wouthuysen transformation, as it has been done in [20] • 
The total energy of the hadron is obtained by adding 
the individual contributions of the jjuarks. The spurious 
center of mass (CM) is corrected [27[ by using intrinsic 
II. THE RELATIVISTIC HARMONIC MODEL operators for the £\ rf and £V V i terms appearing in 

the Hamiltonian. This amounts to just subtracting the 
CM motion zero contribution from the E 2 expression. 

We come now to the description of the quark-antiquark 
potential; it is given by the sum of OGEP and III poten- 
tial 



In the RHM 20], quarks in a hadron are confined 
through the action of a Lorentz scalar plus a vector har- 
monic oscillator potential 



V coni (r) = i(l + 7 o) A 2 r 2 + M 



where 70 is the Dirac matrix 



7o = 
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(1) 



(2) 



Vq( r ij ) = VoGEp(rij) + Vm{nj), 



(10) 



with Tij the inter-quark distance. Among the several 
versions of the OGEP, we have used the following one, 
first derived in 8] from the QCD Lagrangian in the non- 
relativistic limit 



M is the quark mass parameter and A 2 the confinement 
strength. They have a different value for each quark 
flavour. In the RHM, the confined single quark wave 
function ^ is given by 



VbGEp(nj) = -^Aj ■ \j 



(11) 
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with the normalization 
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where the first term is the residual Coulomb energy and 
the second-term the chromo-magnetic interaction leading 
to the hyperfine splittings. Here, is the generator of 
the color SU(3) group for the i th quark, Cj is the Pauli 
spin operator. 

To model the III, we have used the form given in 0, 
HI: 



E is an eigenvalue of the single particle Dirac equation 
with the interaction potential given by (TT|). We can 
perform a unitary transformation to eliminate the lower 
component of 4" such that 
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1 = 1, 

I = 1/2 
S(nj)SsfiS Lt o, 1 = 0. 
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u * = 
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In the above expression, S, L, I are respectively the 
spin, the relative orbital angular momentum and the 
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isospin of the system. The g and g' are dimensioned 
coupling constants. The Dirac delta function appearing 
in ([12]) needs to be regularized for practical calculations. 
As mRef.0,13, 

we have chosen a Gaussian-like function 



The tensor term of III is 
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exp[ 
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The non-central part of the OGEP consists of the 
tensor term ^ogep(^j) ano - the spin-orbit interaction 
Vnnu.pffii)' There are several versions of tensor term in 
literature [26| . We use the expression derived in [8| from 
the QCD lagrangian in the non-relativistic limit and used 
subsequently by many authors (for example [281129^) 



The term Vts(f»j) is responsible for the splitting of the 
3 Lj states with J = L — 1,L,L + 1. The term Vla^j) 
couples states 1 L,/=l and 3 £j=l and due to mass de- 
pendence this term is inoperative when the quarks are 
identical. It is to be noted that the III and the OGEP 
have the same spin dependence except for the Vla term. 



III. FITTING PROCEDURE 
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where Sy = 3((?j ■ f) (Bj ■ f) — o*j ■ Bj . The tensor potential 
is a scalar which is obtained by contracting two second 
rank tensors. Here, f = fi — fj is the unit vector in the 
direction of f. Note that in the presence of the tensor 
interaction, L is no longer a good quantum number. 
The spin-orbit (SO) interaction of the OGEP is given 

by 



The parameters of the RHM are the masses of the 
quarks, M u = Md and M s , the respective confinement 
strengths, A 2 = A 2 ,, A 2 , and the oscillator size param- 
eter b n (= l/0„). They have been chosen to reproduce 
various nucleon's properties: the root mean square charge 
radius, the magnetic moment and the ratio of the axial 
coupling to the vector coupling [20]. The confinement 
strength A Ui d is fixed by the stability condition for the 
nucleon mass against the variation of the size parameter 

K 



OGEPV'yJ — ^ A i A J 



(15) 
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where the angular momentum is defined as usual in terms 
of relative position and the relative momentum P^. 
Unlike the tensor force, the spin-orbit force does not mix 
states of different L, since L 2 commutes with L ■ S, L is 
still a constant of motion, but L z is not. 

The spin-orbit term of III is (see Refs. [3, EH) given by 

vi? (fa ) = V LS (r tj )L -S + V LA [nj )L - A (16) 

The first term in Eq. (|16p is the traditional symmetric 
spin-orbit term proportional to the operator L ■ S. The 
other term is the anti-symmetric spin-orbit term propor- 
tional to L ■ A, where A = \(B\ — B-i). The radial func- 
tions of Eq. (fTB)) are expressed as 




1 



3 I fe=l 



2 I 

\- exp( 



-rfj hi) 



1 



4 



MM, 

h — 



fc=3 



(vkV^) 3 

e xp(-r,y^_ 2 ) 

(??fe-2\/7i : ) 3 



(17) 
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The value of b n for S wave was 0.77 fm [24J and for 
P wave b n is 0.7 fm. It is to be noted that b n is state 
dependent. The parameters associated with the strange 
quark M s and A 2 , have been fitted in order to repro- 
duce the magnetic moments of the strange baryons, ac- 
cording to the procedure described in [3(| ■ The coupling 
constant a, of OGEP is fixed from S wave meson spec- 
troscopy [24[. The value of a s turns out to be 0.2 for P 
wave mesons, which is compatible with the perturbative 
treatment. The parameters of central part of III are g, 
g' and A, the strength and the range of the interaction 
of III, which were fitted to the experimental masses of 7r 
and K mesons [25| . 

Among the non-central parts of the potential, the hy- 
perfine terms of III has 12 additional strength and size 
parameters k and rj (in Eq. [T7l[l"9f respectively. We are 
able to reproduce the light P wave meson masses with all 
r\ and k\ to kq parameters held fixed and by varying only 
the k*[ and Kg parameters. The values of kj and k% are 
listed in Table HU It is to be noted that for each category 
of meson nonet, k-j is held fixed and only the Kg is varied. 
The non-central terms of OGEP are attractive, whereas 
the strengths of the interaction of III (i.e., k) can have 
both positive and negative values [15| . We are led to the 
conclusion that inclusion of III in the formalism is essen- 
tial. This also enables us to bring down the value of a s 
to 0.2. 

The qq wave function for each meson is expressed 
in terms of oscillator wave functions corresponding to 
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TABLE I: Values of parameters used in our model. 



TABLE III: The pseudo- vector meson masses (in MeV). 



b„ 


0.7 fm 


A/ U , d 


315 MeV 


M B 


450 MeV 


a a 


0.2 


Al = Al 


3754.7 MeV fm~ 2 


Al 


3367.45 MeV fm" 2 


9 


0.0847 x 10~ 4 MeV -2 


g' 


0.0535 x 10~ 4 MeV -2 


A 


0.35 fm 


Vi 


0.194 fm 


V2 


0.294 fm 


V3 


0.112 fm 


<H 


u.oui im 


Kl 


-2.213 


K2 


0.191 


K3 


-2.13 


Ki 


2.651 


K5 


20.84 


Kg 


26.43 



TABLE II: Values of K7 and Kg parameters used in our model. 



N 2S+1 Lj 


Meson 




Kg 


l 3 Po 


fa (600) (or a) 


-31.63 


23.92 


l 3 Po 


/o(980) 


-31.63 


17.35 


l 3 Po 


a (980) 


-31.63 


6.81 


l 3 Po 


#o*(1430) 


-31.63 


-21.19 


l 3 Po 


a (1450) 


-31.63 


-87.45 


l 3 Po 


/o(1500) 


-31.63 


-47.26 


1 3 P1 


ai(1260) 


28.8 


1.65 


1 3 P! 


#i(1270) 


28.8 


-1.67 


l 3 Pl 


/i(1285) 


28.8 


-10.32 


l 3 Pl 


/i(1420) 


28.8 


18.59 


1 3 P 2 


/ 2 (1270) 


22.84 


-37.95 


1 3 P 2 


a 2 (1320) 


22.84 


25.52 


1 3 P 2 


# 2 *(1430) 


22.84 


35.19 


1 3 P 2 


/a (1525) 


22.84 


47.23 



the CM and relative co-ordinates. The oscillator quan- 
tum number for the CM wave functions are restricted to 
N cm — 0. The Hilbert space of relative wave functions is 
truncated at radial quantum number n = 2. The Hamil- 
tonian matrix is constructed for each meson separately 
in the basis states of |iV cm = 0, L cm = 0; N 2S+1 Lj) and 
diagonalised. 



Meson 


6i(1235) 


hi (1380) 


#i(1400) 


Experimental Mass 


1229.5±3.2 


1386±19 


1402±7 


Calculated Mass 


1229.22 


1386.42 


1408.34 



TABLE IV: The triplet meson masses (in MeV). 



N 2S+1 Lj 


Meson 


Experimental Mass Calculated Mass 


l 3 Po 


/o(600)(orcr) 


400-1200 


703.54 


l 3 Po 


/o(980) 


980 ± 10 


989.72 


l 3 Po 


a (980) 


984.7±1.2 


986.84 


l 3 Po 


#o(1430) 


1412±6 


1412.07 


l 3 Po 


a (1450) 


1474±19 


1471.77 


l 3 Po 


/o(1500) 


1507±5 


1507.38 


l 3 Pi 


ai(1260) 


1230±40 


1229.54 


l 3 Pi 


#i(1270) 


1273±7 


1272.14 


l 3 Pi 


/i(1285) 


1281.8±0.6 


1281.15 


l 3 Pi 


/i(1420) 


1426.3±0.9 


1426.74 


1 3 P 2 


/a (1270) 


1275. 4±1. 2 


1276.83 


1 3 P 2 


a 2 (1320) 


1318±0.6 


1317.80 


1 3 P 2 


#1(1430) 


1425. 6±1. 5 


1423.61 


1 3 P 2 


/a(1525) 


1525±5 


1525.9 



teraction is not sufficient to reproduce the masses of the 
mesons. If OGEP is taken as the only source of hyper- 
fine interaction, the value of a s necessary to reproduce 
the hadrons spectrum is generally much larger than one; 
this leads to a large spin-orbit interaction, which destroys 
the overall fit to the spectrum. The inclusion of III will 
diminish the relative importance of OGEP for the hyper- 
fine splittings. The important role played by the III in 
reproducing the masses of these mesons (as shown in Ta- 
ble [IV| can be gauged by examining the Table IVl where 
the masses of the scalar mesons calculated after switch- 
ing off the III in the full Hamiltonian are tabulated. This 
is because the tensor and spin-orbit terms of OGEP are 
attractive and hence bring down the masses of the triplet 
state. Also, it is important to note that, the spin-orbit 
terms of III are very weak [15| . The dominant contri- 
bution to the splitting of masses of axial vector mesons 
comes from the tensor term of III, which involves the 
parameters K7 and Kg. It was necessary to tune K7 and 



IV. RESULTS AND DISCUSSION 

The masses of the singlet and triplet P-wave mesons 
after diagonalisation in harmonic oscillator basis with 
n m ax=2 are listed in Table ITnl and llVl respectively. 

Our results indicates that only the OGEP hyperfine in- 



TABLE V: The masses of scalar mesons (in MeV) after diag 
onalisation without III 



Meson 


Experimental mass 


Calculated Mass 


a (1450) 


1474±19 


1002.65 


# *(1430) 


1412 ±6 


1157.45 


/o(1500) 


1507±5 


1259.84 
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Kg, parameters so as to get a reasonably good agreement 
with the experimental masses. Hence, in our model we 
have only two free parameters kj and Kg. Also, the re- 
sults indicate that tensor part of III is crucial and plays a 
dominant role in explaining the masses of P wave mesons 
which is an important result of our investigation. The at- 
tractive or repulsive nature of III being governed by the 
sign of the k. Thus by tuning the k parameters appro- 
priately, we are able to reproduce the meson masses in 
our model. 



A. Pseudo-vector meson nonet (I 1 Pi) 

We have investigated three mesons of the l*Pi pseu- 
dovector meson nonet with J = 1^ , namely, 
6i(1235), hi (1380), # x (1400) [p. It may be pointed 
out here that there is no contribution from the III for 
the singlet states except for the 1 P\ state in the K-sector. 
In the K-sector, the singlet P state receives a significant 
repulsive contribution of 129.75 MeV from the off diago- 
nal matrix element (ME) of ( 3 -Pi|Vla| 1 Pi)- The masses 
calculated by our model are tabulated against the corre- 
sponding masses [3l| in Table IIIII 

B. Scalar meson nonet(l 3 P ) 

The spectrum of the scalar meson nonet is very large 
and the actual number of resonances in the region of 1-2 
GeV far exceeds the number of states which the conven- 
tional quark models can accommodate. Several of these 
states, however, have been interpreted as exotic mesons. 
It is well known that a qq meson with orbital angular mo- 
mentum I and total spin s must have parity P = (— l) i+1 
and charge conjugation quantum number C = (— l) l+s . 
On this basis, we define an exotic meson to be one which 
does not have the above spectroscopic configurations. 
Thus a resonance with J PC = , + ~, 1 _+ , 2+~ ■ ■ • 
are exotic. Such states could be a gluonic excitation such 
as a hybrid (qqg) or glue ball (2g, 3g ■ ■ ■ ) or a multi quark 
staXc{qqqq) . 

The particle data group(PDG) lists isoscalar states, 
the a (980) and a (1450) [Si having masses of 984. 7 ±1.2 
MeV and 1450 ±40 MeV respectively. Theories based on 
chiral sigma models with three flavors [HJ suggest that 
ao(980) would form a scalar nonet. The scalar ^(1430) 
is well established. Several groups have claimed dif- 
ferent isoscalar structures close to 1500 MeV [33|, 1341 ]. 
In this work, we focus our attention only on the non- 
exotic scalar mesons with J PC = ++ and assigned 
as jb(600V_/o(980), a (980), 2^(1430), a (1450) and 
/o(1500) [3l|.The isosinglet /o(980)and and isotriplet 
ao (980) lower than 1 GeV are considered as a chiral part- 
ner (cr nonet) of the pseudo scalar ((it nonet) in SU(3) 
chiral symmetry. The light scalar mesons with masses less 
than 1 GeV are not considered as a conventional L=l 
scalar nonet by many authors. 



C. Axial vector meson nonet (l 3 Pi) 

In our model, for axial vector mesons, the tensor and 
L ■ S parts of OGEP and III have opposite signs. The 
contributions due to tensor terms are repulsive, whereas 
those due to L-S terms are attractive. As the OGEP has 
the same strength parameter for these terms, the contri- 
bution of the hyperfine interaction terms of OGEP is neg- 
ligible whereas, due to the different strength parameters 
Ki, the corresponding terms of III contribute differently. 
Besides, the contribution of III to the masses is also sig- 
nificant because of the different radial form of tensor and 
spin-orbit terms. We have treated free parameters 
so as to reproduce the masses of ai(1260), ifi(1270), 
/i(1285) and /i(1420). However, it should be noted 
that the ai(1260), with 1=1 has a significant width of 
400 MeV and has a dominant decay channel ai — > pir. 
This property makes the determination of its mass dif- 
ficult. The QCD sum rules [3|| produce a mass of 
1150±40 MeV. According to Bowler [36|], the a\ mass and 
width are safely within the ranges ~ 1235 ± 40 MeV and 
400 ± 100 MeV respectively. These values are in agree- 
ment with those currently adopted by PDG [3l[, i.e., 
mass of 1230 ± 40MeV and width 250 MeV to 600 MeV. 
In the K-sector, we have fitted to iV(1270). The contri- 
bution from the matrix element (ME) (1 1 Pi|Vla|1 3 Pi) 
has been found to be significant. PDG cite two /i me- 
son states [3l| with J PC = 1 ++ , namely, /i(1285) and 
/i(1420). There has been considerable discussion on the 
quark structure of these mesons (3?| • We have been able 
to fit the masses reasonably well as shown in Table IIVI 



D. Tensor meson nonet (1 3 P 2 ) 

We consider some of the well established members of 
the tensor meson nonets, with J PC = 2 ++ , i.e., /2(1270), 
a 2 (1320), #|(1430) and /£(1525). The contributions due 
to tensor and L ■ S terms of OGEP and III bear opposite 
signs. The tensor potential is attractive whereas L ■ S 
part is repulsive. However, the off diagonal tensor ME 
(3P2|V ( J GE p| 3 F2) is strongly repulsive. In our model the 
mass difference between fi and / 2 essentially comes from 
the off diagonal ME of tensor potential of OGEP and III. 

In literature some more 

jpc _ 2 ++ states like 
/ 2 (1520), / 2 (1810), / 2 (2010), / 2 (2340) have been con- 
sidered. Of these, / 2 (1810) is likely to be the 2 3 P 2 state 
[13 ■ Our model prediction for / 2 (1810) is 1724.52 MeV. 



V. CONCLUSIONS 

The mass spectrum of P wave mesons is considered in 
the frame work of RHM with the conventional OGEP 
and by including III. The inclusion of III consequently 
diminish the relative importance of OGEP. The III also 
restricts a s to be 0.2 and thus justifying the perturba- 
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tive truncation of multi gluon exchanges. The near mass 
degeneracy of the experimentally established iso-doublct 
states of the scalar and tensor meson nonets Kq and 
K2 could be accounted by the off diagonal tensor ME 
of OGEP and III. The simultaneous mass degeneracy 
of the pseudo-vector K\b and axial vector K\a which 
mix to give physical i^i (1270) and if 1 (1400) states ob- 
served experimentally could be accounted for by the anti- 
symmetric spin-orbit term Vla of III. As we have shown, 
RHM with OGEP and III provides a quite good descrip- 
tion of the pseudo-vector, scalar, axial vector and tensor 
P-wave mesons with the same constituent quark masses, 
oscillator size and OGEP strength a s . The calculated 
masses of the P wave mesons are in good agreement with 



the experimental masses. 
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